Spin-orbit locking as a protection mechanism of the odd-parity superconducting state 

against disorder 
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Unconventional superconductors host a plethora of interesting physical phenomena. However, 
the standard theory of superconductivity suggests that unconventional pairing is highly sensitive to 
disorder, and hence can only be observed in ultraclean systems. We find that due to an emergent 
chiral symmetry, spin-orbital locking can parametrically suppress pair decoherence introduced by 
impurity scattering in odd-parity superconductors. Our work demonstrates that disorder is not an 
obstacle to realize odd-parity superconductivity in materials with strong spin-orbit coupling. 

PACS numbers: 74.20.Rp, 73.43.-f, 74.20.Mn, 74.45. -|-c 



The quest for topological superconductors [TWSl is an 
exciting research area in condensed matter physics. A 
necessary condition for topological superconductors is 
unconventional pairing symmetry. Based on an ear- 
lier parity criterion [6], it has been shown that under 
fairly general conditions, time-reversal-invariant topolog- 
ical superconductivity is realized if the pairing symme- 
try is odd under spatial inversion [21 [5] . Such odd-parity 
pairing occurs in the Ballian-Werthamer phase of super- 
fluid helium-3 ^, and likely in certain heavy fermion 
superconductors p^OlITT]. 

A recent theoretical study 2\ suggests that doped 
narrow-gap semiconductors are candidates for odd-parity 
topological superconductors. Here the strong spin-orbital 
coupling in the band structure favors a novel interorbital, 
odd-parity pairing, even when the mechanism for super- 
conductivity is conventional clcctron-phonon interaction. 
Experimentally, several materials in this class, including 
Cu-doped BizScg [T^, Tl-doped PbTe [TS], and In-doped 
SnTe ;14j , exhibit superconductivity with unusually high 
transition temperatures (2 — 4K) relative to their low car- 
rier densities (~ 10^*'cm~'^). Recently, some evidence of 
non s-wave pairing in Cua;Bi2Se3 has been reported [TSl- 
[T7] . including the presence of a zero-bias conductance 
peak in point-contact spectroscopy [T8j [19] . 

An important issue in the study of unconventional 
superconductors is their robustness against disorder. 
Within the Bardeen-Cooper-Schrieffer theory, s-wave 
pairing is immune to nonmagnetic impurities [20], while 
other pairing symmetries are more fragile [21]. For in- 
stance, the transition temperature of Sr2Ru04 that is 
believed to be a spin-triplet superconductor is strongly 
suppressed by disorder[32]. This seems to be a major ob- 
stacle for realizing odd-parity superconductivity in doped 
narrow-gap semiconductors. 

In this Letter, we study the effect of disorder on the 
proposed odd-parity superconducting state in narrow- 
gap semiconductors [H |2]. We mostly focus on scalar 
impurities which are usually the most common type of 
disorder. Such an isotropic scattering potential arises 
from Coulomb interactions between the electrons and lat- 



tice defects. Contrary to conventional wisdom, we find 
that the destructive pair-breaking effect of disorder is 
dramatically suppressed by an approximate chiral sym- 
metry in the spin-orbital locked band structure. In view 
of our study, the prospect of topological superconductiv- 
ity in narrow-gap semiconductors appears brighter than 
before. 

For concreteness, we start our analysis with the follow- 
ing non-interacting four band k • p Hamiltonian: 



-f^o(k)= ^"^(k) [mn,j.+vayik,^Sy 



-kySx) + v^k^ay\tl;{k) 



(1) 



Hq has the form of a massive Dirac Hamiltonian. Here, v 
and Vz are the velocities of the electrons in the xy plane 
and z direction, respectively, while the Dirac mass m 
(not to be confused with the effective mass) determines 
the energy gap between the bands. Eq. (II]) describes 
the band structure of a broad class of narrow-gap semi- 
conductors with inversion symmetry near time-reversal- 
invariant momenta [5]. For example, in the context of 
doped Bi2Se3, PbTe and SnTe, Si are Pauli matrices in 
spin space, and Ui are Pauli matrices associated with the 
orbital degrees of freedom [ij [J^ . Specifically for doped 
Bi2Se3, we use cr^ = ±1 to label the two p^-like orbitals 
located at the upper and lower part of the quintuple layer 
unit cell. The two orbitals (+ and — ) transform into each 
other under spatial inversion with respect to the center 
of the unit cell. 

Short range (phonon mediated) attractive interactions 
can generate pairing with two distinct symmetries. One 
is the conventional spin-singlet intraorbital state, while 
the other is the unconventional spin-triplet, orbital sin- 
glet paired state [2]: 

Ai cx^(V'.,t(k)V.4(-k) - ^,,;(k)^,,^(-k)); (2a) 

cr=± 



A2 cx^a(^.,^(k)^_„,;(-k) + ^.,;(k)V'-.,t(-k)). (2b) 
<j=± 

The Ai pairing is invariant under all symmetry opera- 
tions of the D3d point group (s-wave) , whereas A2 is odd 



under spatial inversion (which interchanges the two or- 
bitals) and belongs to the Aiu representation. In Ref. 2 
it has been shown that, in the absence of disorder, the 
strength of the interorbital and intraorbital attractive in- 
teraction determines which pairing susceptibility diverges 
more strongly. Here, we examine the effect of disorder 
on the odd-parity pairing assuming the interaction fa- 
vors this state. In contrast to the well known p-wave 
spin-triplet superconductivity, here the A2 order param- 
eter is independent of momentum, but has a nontrivial 
internal structure in orbital and spin space. In the tradi- 
tional p-wave superconductors the order parameter varies 
over the Fermi surface, and scattering of a Cooper pair 
from the states k, — k into k', — k' results in their deco- 
herence. Consequently, superconductivity dies when the 
elastic scattering rate becomes comparable to the order 
parameter, 1/r — >■ A. Our result shows that decoherence 
effects in the spin-triplet, orbital-singlet paired state are 
significantly suppressed, and superconductivity survives 
much stronger levels of disorder until {m/ ^Y jr — > A. 

To show quantitatively the effect of disorder on the 
onset of superconductivity, we examine the pairing sus- 
ceptibility near the transition. First, we transform into 
the eigenstate basis, in which the Hamiltonian is diag- 



onal. 



^o=Ek,,^(k) 



where 
w?A;?. De- 



cT(k)c,(k)-d](k)d,(k) 

the dispersion is i?(k) = ., ra^ + v'^{k'^ + kT) 
spite the strong spin-orbital mixing in Hq, the presence 
of both time-reversal and inversion symmetries protects 
the two-fold degeneracy of the upper (ci, C2) and lower 
(di, (^2) bands at every k. We chose to work in a par- 
ticular basis [^, which we dub "pseudo-chiral", where 



^v(k,C^n) = ^'5= 



the states at k and — k are related (up to a k-dependent 
phase factor) by time reversal (0) and inversion [P) op- 
erations: 

eci(k)e-i - ci(-k), ec2(k)e-i - c2(-k) 

Pci(k)p-1 ^ C2(-k), Pc2(k)p-1 ^ Ci(-k). (3) 

As long as the Hamiltonian has both time-reversal and 
inversion symmetries, it is always possible to use the ba- 
sis satisfying the transformation properties described in 
Eq.[3l 

In the pseudo-chiral basis, the superconducting order 
parameters (2a) and (2b) are given by: 



^k',k y-kM' 
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k',k 
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FIG. 1. The effect of disorder on the single particle Green's 
function and pair susceptibility. The seff-consistent equation 
for the self-energy (Eq. u\ is illustrated in (a) . The double line 
represents the dressed electron Green's function, and the dashed 
line denotes the impurities. The equation for the Cooperon 
(Eq. |9| is illustrated in (b). Here, the Cooperon is drawn as a 
filled triangle vertex. 






Ai(k)(xe-^* -^ {-ly [(c,(k)c,(-k)) + (d,(k)d,(-k))] 



(4a) 



A2 (k) (X e-^^ cos ak^ [(c, (k)c, (-k)) - (d, (k)d, (-k))] 

- 2e-'*sinak [(ci(k)d2(~k)) + (c2(k)di(-k))] . (4b) 

In the derivation of the above expressions and from now 
on we rescaled the z component of the momentum by its 
velocities, k = {kx, ky, Vzkz/v). While the expressions for 
the order parameters contain the azimuthal angle be- 
tween (kx, ky) and the x axis, they are independent of the 
polar angle 9 between k and the z-axis. The parameter 
ak — sin^ {m/ECk)) is a consequence of the interorbital 
mixing by the max term. 

In the limit m = (ak = 0), the k • p Hamiltonian has 
a U{1) chiral symmetry, [Hq, o-ySz] — 0. Thus, the energy 
eigenstates labeled by j = 1(2) have a well-defined chi- 
rality -|-1(— 1), which is evident from Eq. HI Importantly, 
pairing in both the even- and odd-parity states only oc- 
cur between electrons of the same chirality, but the order 
parameters differ by the relative phase between Cooper 
pairs of opposite chirality (of different label j), tt for Ai 
and for A2(2]. Since scalar disorder potential (which 
is insensitive to orbitals and spins) does not break chi- 
ral symmetry, impurities can scatter electrons only be- 
tween states of equal chirality. Hence, the disorder af- 
fects the even-parity pairing Ai and odd-parity pairing 
A2 in the same way. Then, it follows from the Anderson 
theorem [50] that as long as the system is far from local- 
ization, both pairings are completely robust against dis- 
order. The magnitudes of two order parameters at zero 
temperature differ only by the strength of the pairing in- 
teraction in the two channels, A^^i 2 = r2exp{— l/fpAf} 
with z/Q the density of states at the Fermi energy, Xi the 
attractive interaction in the even- or odd-parity channel, 
and rj an ultraviolet cutoff. 

For nonzero m in the Hamiltonian Hq, the Bloch wave- 
functions are no longer chiral eigenstates, and the similar- 
ity between the two order parameters is broken. Without 
loss of generality, we assume that the chemical potential 



lies in the upper energy band, ^ > \m\. Since only elec- 
trons on the Fermi surface contribute to the divergent 
pairing susceptibility at low temperature, the effective 
pairing order parameter involves only the c electrons of 
the conduction bands: 

Ai(k) ex e'^ [(ci(k)ci(-k)) - (C2(k)c2(-k))] (5a) 



A2(k) ex e'-^cosa [(ci(k)ci(-k)) + (c2(k)c2(-k))] . (5b) 

Here a = sin~ (m/^) is evaluated at the Fermi en- 
ergy. One can see that for m ^ the effective pair- 
ing potential in the odd-parity channel between elec- 
trons in the conduction bands is reduced due to mix- 
Assuming 



equation: 



ing with the valence electrons (see Eq. 4b I. 
the strength of the inter-orbital attraction is indepen- 
dent of the parameter m//^, the magnitude of the odd- 
pairing order parameter at zero temperature becomes 
A2 = 1^ exp{ — 1/z^Q cos^ aXi}. The order parameter A2 is 
reduced from its ?n = value not only due to the depen- 
dence of the density of states on m//z, Va = vq cos a, but 
mainly because the attractive interaction is weaker by a 
factor of cos a. We wish to emphasize that this change 
in the odd-symmetry order parameter is not caused by 
pair breaking. Pair breaking effects reduce the transi- 
tion temperature, but not the zero temperature order 
parameter |26j. 

Introducing scattering by a scalar disorder poten- 
tial adds a term into the Hamiltonian which is non- 
diagonal in the momentum, _ff(k, k') = iJo(k)(5k,k' + 
VinipV'^(k)^(k). Transforming into the pseudo-chiral ba- 
sis, one can see that the impurities mix all bands. How- 
ever, in the limit epT 3> 1, all leading order processes 
occur near the Fermi surface, and we can again restrict 
our attention to the conduction bands. The matrix ele- 
ment between two states at the Fermi energy, il^i{k) and 
'(/'j(k'), is given by 



, _ , ^(k,k') A(k, 
''k.k' - I 4(-k,k')sina A{- 



-k') 
k,- 



sma 
■k') 



(6) 



|e*("^-'^')cosf cosf^ 



Here ^(k,k') = ^i^p^ 
is equal to the wavefunction overlap between two spins 
pointing along k and k' directions. In the limit m = Q 
[a — 0) the off-diagonal matrix elements vanish, restor- 
ing the chiral symmetry 

To find the corrections to the electron self-energy S 
due to scattering by disorder we use the self-consistent 
Born approximation. Then the self-energy matrix for the 
two conduction bands is diagonal and determined by two 
processes, intra-band scattering and scattering between 



the two conduction bands (for illustration see Fig.|l(a) I 



■^k,Lj„ 



I 

dk' 



sgn(w) 



dk' 

W) 



3 I '^k.k' 



'■G 



k' LJ„ 



(2^) 



iJ^'iyi^'i.iQJ 



3 Mc,k' ""k'.k 



(7) 



where G^ ^(k,a;„) 



lUJr, 



{E{k) - y) - I](k,w„) is 



the Matsubara Green's function of the electrons in the 
band i. Examining the above equation, one can see 
that the two scattering processes interfere constructively. 
Consequently, the elastic scattering time satisfies 1 = 

-KVoiV^^pTil + SIT? a). 

As we explained above, while the transition tempera- 
ture is highly sensitive to pair decoherence mechanisms, 
the zero temperature order parameter is not modified. 
Therefore, we calculate the transition temperature into 
the superconducting phase for the even- and odd-pairing 
states from the corresponding pairing susceptibilities of 
the normal state: 



/ = AT ^ Tlpi, 



G-' 



Cij-,k,m{(^n)^k,m- (8) 



Here, Tfj — Si,j[Si^i — 6.1^2] for the spin-singlet pairing 
i = 1, and rf ■ = Sij cos a[Si I + J^^] for the orbital- 
singlet pairing £ = 2. The matrix Ci,j-k,mi'-^n) is the 
Cooperon describing multiple scattering events of two 
electrons (a Cooper pair) in the particle-particle channel. 
A pole in the Cooperon at (u;„ = means that scattering 
by impurities does not result in decoherence of Cooper 
pairs. In other words, the probability of an electron in 
state k to be scattered into state k' is equal to the prob- 
ability of its partner to be scattered into the partner of 
k'. A Cooperon with a finite mass, on the other hand, 
implies that pairing is suppressed by disorder and that in 
the superconducting state there are sub-gap excitations. 
The Cooperon can be expressed in terms of the single 
particle Green's function and impurity scattering poten- 
tial: 



dk 



^ J (27r) 



3 '^k,k''^-k.-k'<^k' 



(9) 



The self energy is found from the following self-consistent 



The only four components of the Cooperon that enter 
the susceptibilities are Ci^i-jj. In the absence of impuri- 
ties Ci^i;2.2 — (^2.2:1,1 = and only the two components 
Ci,i;i,i = C'2,2;2,2 matter. As a result, there is no differ- 
ence in the effect of disorder on both order parameters. 
Similar to the single-particle self-energy, the Cooperon 
includes processes in which the electrons remain in the 
same band after scattering and those in which at least 
one electron changes its band. However, while the inter- 
ference between the scattering events that determine the 
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transition temperature enters the Cooperon as a mass 
term eliminating its divergence. When 1/t^ is larger 
than the order parameter in the absence of pair deco- 
herence, superconductivity disappears. We found that 
disorder introduces a pair breaking mechanism for the 
odd-parity pairing but not for the even-parity pairing, 
and that the only difference between the two order pa- 
rameters is the phase between the condensates in the two 
conducting bands. Therefore, we conclude that the deco- 
herence is due to the tendency of the impurity scattering 
to favor a relative phase of tt between the conduction 
bands (the singlet pairing), suppressing the spin-triplet 
transition temperature: 



In 



TO ^ I 2 



"^ ' 2 ^ 2ttT,{t)t^ 



(12) 



FIG. 2. The critical temperature as a function of the level of 
disorder TvugVi^p for various values of a. Both Tc and Trj^oV^mp 
are given in units of of T^, the critical temperature for Fimp = 0. 
For comparison, the dashed line shows the transition tempera- 
ture into the odd-parity superconducting state in the absence 
of spin-orbital locking (as in Helium-III). One can see that up 
to a ~ 0.5 the critical level of disorder in which the unconven- 
tional order parameter disappears is dramatically higher than 
in Helium-III. At a < 0.3 the effect of moderate disorder on T^ 
is almost unnoticeable. 



elastic scattering time is constructive, the Cooperon is a 
sum of constructive and destructive interference terms. 
This can be seen in the expressions for Ci^i^jj: 



a 



l + 2|w„|r 
-I- (-1)'+^' 



1 + 2\uj„\t 



imp 

1 



sin^ a) 



(10) 



2\uj„ 



1 + 2|a;„|T - T:VaV^^pT{l -t- sin^ a 



Note that the above expression for the Cooperon has 
been calculated only assuming that epT ^ 1. Now the 
two susceptibilities, x^ and y^ , are no longer identical. 
The even-pairing susceptibility includes only the con- 
structive interference processes, x^ = 47rr^^ i^aT[l — 
7rr.„y4pr(l + sin2a)(l-2|a;„|r)]-i = 2^TE J^a/|^„|, 
and is clearly insensitive to disorder. In contrast, the 
susceptibility for odd-pairing is determined by the de- 
structive interference, and hence, affected by impurities: 



X 



T 



•^—^ irVa cos'^ a 

^-^ \ijJn\ + sin^ aJTiX -f sin^ a) 



(11) 



21^0 (cos a) 



In- 



27rT 



^ 



. 2 

sm a 



2TrT{l 



sm a}T 



where ipix) is the digamma function. 

Let us explain the above result. The pair decoherence 
rate, l/r^p = (sin a)/{l + sin a)T, which reduces the 



Here we use T^ to denote the transition temperature at 
1/r = and finite a, and C — -k I2e^ with 7 the Euler 
constant. The suppression of T^ as a function of disorder 
for various values of a is illustrated in Fig. [2] 

To understand better the peculiarity of this result, 
it is instructive to compare our result, applicable for 
the narrow-gap semiconductors, with the BW super- 
fluid phase in Hc-III. Although the latter has a single 
band without spin-orbit coupling, one can choose a basis 
in which spin is locked to be parallel, (/>i(k), or anti- 
parallel, (/'2(k) to the momentum. This basis resem- 
bles the chiral basis we used here. Correspondingly, in 
both cases the odd-parity order parameters can be writ- 
ten as e^'^[(/)i(k)(/.i(-k) -I- 02(k)(/)2(-k)]. The key dif- 
ference between systems with and without spin-orbital 
locking manifests itself in the impurity scattering. While 
in Helium-III , the matrix element for impurity scatter- 
ing between 0i(k) and (/'2(k') are larger than the diag- 
onal terms, in the problem studied here the intcr-band 
scattering is parametrically smaller than the intra-band 
one by a factor of sin a — m/ ^, due to the approximate 
chiral symmetry that becomes exact as m — )■ 0. Thus, 
the pair decoherence in the narrow-gap semiconductors 
is significantly weaker than in Helium-III. Note that this 
observation is only correct for scalar disorder, and does 
not hold for other types of scattering potentials, such as 
magnetic impurities (ex s) or orbital dependent poten- 
tials ( ex (?). These non-scalar disorder potentials cause 
stronger pair decoherence and suppresses the unconven- 
tional superconducting state. 

Our result can be generalized to 2D bilayer band struc- 
tures obtained by setting fc^ = in the k ■ p Hamilto- 
nian ([I]) . Recently, it has been proposed [3] that this class 
of bilayer systems with Rashba spin-orbit coupling favors 
odd-parity superconductivity similar to the 3D case. Our 
analysis shows a similar robustness against disorder due 
to the chiral symmetry (See Supplemental Material |28|). 

To conclude, we showed that certain types of odd- 
parity pairing in doped narrow-gap semiconductors can 



survive from a fairly large amount of impurity scatter- 
ing. The relative robustness of these systems results 
from an approximate chiral symmetry in the spin-orbital 
locked band structure and the odd-parity pairing order 
parameter. Although scattering by disorder reduces the 
phase coherence of the Cooper pairs, the dephasing rate 
l/r<^ = TTfaVj^pSin^ a vanishes as the Dirac mass {= 
band gap) in the band structure goes to zero, when the 
chiral symmetry becomes exact. Finally, we note that 
in addition to centrosymmetric materials studied in this 
work, strong spin-orbit-coupling in asymmetric interface 
structures can also protect unconventional superconduc- 
tivity against disorder |29j . 
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SUPPLEMENTARY MATERIALS 

Choice of basis 

In the supplementary material part, we elaborate on the 
transformation connecting between the basis of spin (s =t 
, \.) and orbital {a — ±) eigenstates and the pseudo-chiral 
basis. The creation operators of electrons in the conduc- 
tion bands, c-, and valence bands, d^ can be written as: 



c](k) 



rfj(k) 



a,s 

^t;k,,(a,s)<,(k). (13) 



The functions Uk,_, and Vj are four components vectors in 
the spin and orbital space, (+ t, + i, — ti ~ 4-)- For conve- 
nience, we rescale k^, by its velocity, k = {kx,ky,Vzkz/v). 
The transformation functions ui^j to the conduction bands 
are given by 

1 /• 9 -ia/2 ■ S ilrh + a/2) 

Mk,i = — =(icos-e ,sm-e^^ ', 

^ ia/2 ■ ■ 9 i(S-a/2)\ 

COS -e ' ,jsin-e^ ' '), 

1 / . ^ -ia/2 • ^ -i{<i-a,/2) 

Wk,2 = — (sm-e ,icos-e ^^ \ 

■ ■ S ial2 9 _i(rf, + Q/2)^ ., .N 

— jsm-e ' , — cos-e ^^ ' '). (14) 

where 9 and (j> are the polar and azimuth angle of k, re- 
spectively, and a = sm~^ {m/fj,). Similarly, for the valence 
bands: 

1 / 9 -ia/2 ■ ■ 9 ild,+al2) 

Wk,i = — (-cos-e ,-ism-e^'^ ', 

9 ia/2 ■ 9 i(rf,_a/2)N 

— jcos-e ' , — sm-e^ ' '), 



J- / . • " -ia/2 " 

Wk,2 = — =(ism-e , cos-e 

. 9 i^i2 . 9 _ 
— sin -e ' , —I cos -c 



ia/2 O -i((h-a/2) 

i(cl>+a/2)^ /-j^gx 



Using the above transformation we obtained the expres- 
sion for the Hamiltonian, order parameters and the im- 
purity scattering matrix in the pseudo-chiral basis. This 
transformation allowed us to focus on pairing, as well as 
elastic scattering events in the vicinity of the Fermi sur- 
face. 



3D Hamiltonian of Eq. (1) in the main text with k^ — 0: 
-ff2D(k) = V^(k) [max + va^{kxSy — fc^Si,)] ?/)(k). (16) 

Here ctz = ±1 is the layer index and Sz = ±1 is the 
electron spin along the direction perpendicular to the 
planes. As explained in Ref. [Ij, the top layer is subject 
to an electric field created by the bottom layer, and vice 
versa. These fields lead to an intra-layer Rashba spin- 
orbit coupling term with an opposite sign for the two 
layers az{kxSy — kySx)- In the 2D Hamiltonian, the mass 
term generates inter-layer hopping. The above Hamilto- 
nian applies to a wide class of materials including thin 
films/quantum wells of narrow-gap semiconductors like 
Bi2Se3, as well as oxide heterostructures[3J. 
The possible pairing symmetries coincide with those sug- 
gested in Ref. 2 for the 3D Hamiltonian 3 . In particu- 
lar, for fc^ = the orbital-singlet spin-triplet pairing (A2) 
studied in the main text, corresponds to inter-layer spin- 
triplet pairing in the bilayer system. Due to the same chiral 
symmetry discussed in the main text, there is no difference 
in the effect of disorder on A2 in the 2D and 3D cases. 
In addition to the two pairing symmetries discussed in the 
main text, Ref. |2j proposed the following unconventional 
spin-singlet state: 



A3 oc^ cr{^,,t(k)^,4(-k) - V<T,Uk)V'<.,t(-k)); 



(17) 



Unlike the case of A2, the effect of disorder on this paired 
state dramatically changes when dimensionality is reduced 
from three to two dimensions. In the 3D case, A3 has 
nodes, and the paired state is not protected by the chi- 
ral symmetry even when m — 0. As a result, pairing is 
strongly suppressed by disorder. In contrast, in 2D the or- 
der parameter, which corresponds to spin-singlet 2D topo- 
logical superconductor 3_ with opposite phase on the two 
layers, is fully gapped. Moreover, in the limit tti — >■ 0, 
a new chiral symmetry emerges and the Hamiltonian is 
invariant under the transformation ilja;s — >■ e"^ '4'tT,s- For 
9 = 7r/4, the odd-parity order parameter A3 transforms 
into Ai. Thus, as long as this symmetry exists, A3 is pro- 
tected against the effect of both scalar and a^ type of 
impurities. For a finite mass term the suppression of the 
transition temperature into A3 is parameterically weaker 
than it is in the 3D case. It will be interesting to find such 
a s± odd-parity topological superconductor. 



2D bilayer Rashba systems 

In this section we discuss the effect of disorder on un- 
conventional pairing in 2D bilayer systems with Rashba 
spin-orbit coupling. Such a system is described by the the 
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